Entropy of polydisperse chains: solution on the Bethe lattice 
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We consider the entropy of polydisperse chains placed on a lattice. In particular, we study a 
model for equilibrium polymerization, where the polydispersivity is determined by two activities, 
for internal and endpoint monomers of a chain. We solve the problem exactly on a Bethe lattice 
with arbitrary coordination number, obtaining an expression for the entropy as a function of the 
density of monomers and mean molecular weight of the chains. We compare this entropy with the 
one for the monodisperse case, and find that the excess of entropy due to polydispersivity is identical 
to the one obtained for the one-dimensional case. Finally, we obtain an exponential distribution of 
molecular weights. 

PACS numbers: 65.50.+m,05.20.-y 



I. INTRODUCTION 

The problem of the thermodynamic behavior of di- 
atomic molecules which are adsorbed on two-dimensional 
surfaces has attracted the interest of researchers already 
in the thirties In the simplest model of this kind, 
only excluded volume interactions are taken into account, 
and the relevant quantity is the entropy of placing dimers 
(which occupy two first neighboring sites) on the lattice, 
which is a fundamental equation of the system. For the 
particular case of full occupancy and two-dimensional lat- 
tices, this problem was solved exactly by Kasteleyn and 
Temperley and Fisher Q • The model may be generalized 
in many ways. For example, an energy may be associated 
to the configurations of dimers on the lattice, and such 
models may display rather peculiar phase transitions Q . 
Also, we may generalize the athermal model allowing for 
chains of more than two monomers (M-mers) and also 
considering only a fraction p of the sites of the lattice 
occupied by monomers. 

The entropy of Af-mers (chains of molecular weight 
M), as a function of the fraction of lattice sites occupied 
by monomers p, is a fundamental equation of this system, 
and thus contains all thermodynamic information. The 
entropy may be defined as 



sm(p) 



- lim lnT(N p ,M;V), 

V V — >oo 



(i) 



where T(N P , M; V) is the number of ways to place N p 
chains with m monomers in each on the lattice with V 
sites, and the thermodynamic limit is taken with fixed 
density of occupied sites p — N p M/V. Besides the ex- 
act result mentioned above for S2(l) in two dimensions, 



in the literature we may find series expansions estimates 
[3], closed form approximations (Bethe and Husimi lat- 
tices) and transfer matrix calculations with finite size 
extrapolations [6[ on the square lattice. 

It is frequent in polymeric systems that the chains do 
not all have the same number of monomers, that is, a 
polydisperse set of chains is found. One particular model 
for polydisperse chains is the equilibrium polymeriza- 
tion model proposed by Wheeler, Kennedy and Pfeuty 
0, [H and applied to the study of equilibrium polymer- 
ization of sulfur. The system is is defined on a lattice 
(celular model), so that each of the V cells is occupied 
by a monomer (a Ss ring in the case of sulfur). Each 
monomer may be active or inactive (open and closed 
rings for sulfur), and active monomers in adjacent sites 
may connect, forming linear polymers which are self- and 
mutually avoiding walks. Active monomers which are 
not connected to any other ones are considered one-site 
polymers. The statistical weight of chain composed by 
m monomers (and therefore by m — 1 bonds) is given by 



A'i if M = 1, 
2K 1 (K' ) M - 1 if M > 1. 



(2) 
(3) 
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The additional factor 2 in the weight of multi-site chains 
is necessary for the exact correspondence of the equilib- 
rium polymerization model to the n-vector of magnetism 
in the formal limit n — > 0, but it may also be justified 
by a combinatorial argument 0. This model, as well as 
related ones, was studied in some detail in the literature 
[9(, and experimental realizations of this situation are 
found in equilibrium polymerization of sulfur and in the 
so called living polymers [l(| [H[ ■ It is expected that for 
the polydisperse case the entropy is higher as compared 
to the monodisperse case, due to the additional degrees 
of freedom. Another relevant point is the distribution of 
molecular weights M. Recently the entropy of polydis- 
perse chains and the distribution of molecular weights of 
the chains were calculated in the one-dimensional case 



2 



|12j |. In this paper we address these questions for the 
solution of the model on the Bethe lattice with coordi- 
nation number q. As in the one-dimensional case, on 
the Bethe lattice it is also possible to reach closed form 
expressions for the entropy and for the distribution of 
molecular weights. 

I the section Ql] the problem is defined in more detail 
and the relation of the equilibrium polymerization model 
to a simpler version, without the presence of one-site 
polymers, is shown. The entropy and the distribution 
of molecular weights on the Bethe lattice is obtained. 
Final discussions and comments may be found in section 

En 



II. DEFINITION OF THE MODEL AND 
SOLUTION ON THE BETHE LATTICE 



where the sum over configurations is now restricted to 
chains with more than one monomers, the variable x\ is 
equal to 1/2 and the activities of endpoint and internal 
monomers are 



and 



1 + 2K lXl 



1 + 2K x x x 



(7) 



(8) 



respectively. We may now restrict our discussion to 
the model without one-site monomers, described by the 
grand-canonical partition function E(z e , zf, V). In Fig. [1] 
a possible configuration of the chains is shown. If the lat- 



The partition function for model for equilibrium poly- 
merization described above may be written as: 



Y(K X ,K P , x V ,V) = ^2(2Ki) N ' (K' p ) 



(4) 



where the sum is over all possible configurations of the 
polymers on the lattice, N p is the number of polymers 
(including one-site chains) of the configuration, Nb cor- 
responds to the number of bonds and X\ = 1/2 is a fac- 
tor which assures the proper counting of configurations, 
as discussed above. We may now perform a partial sum 
over configurations in this partition function. Let us con- 
sider a particular configuration of the polymers which 
occupy more than one site (M > 1). The number of 
monomers in internal site of the chains will be called Ni 
and the number of the monomers on endpoints of the 
chains is N e . The number of bonds in this configura- 
tion is Nf, = Ni + N e /2 the total number of chains is 
N p = N% + N e /2 and the number of chains with more 
than one monomer is equal to N e /2. We may now write 
the partition function as: 



Y{K u K p ,xi;V) = 
^'(2if 1 )^ /2 K) JVe/2+Ari E"(2X 1 a ;i ) Wl : 



(5) 



where the first sum is restricted to configurations of 
multi-site chains and the second sum is over the configu- 
rations of the one-site polymers for a fixed configuration 
of the other chains. It is now trivial to perform the sec- 
ond sum, which is equal to (1 + 2K\X\) V ~ Ni ~ N " , where 
we notice that V — Ni — N e is the number of sites not 
occupied by monomers in multi-site chains and each of 
these sites may be either empty (weight 1) or occupied by 
a one-site polymer (weight 2K\ X \). Thus, the partition 
function of the equilibrium polymerization model may be 
written as: 

Y{K x ,K p ,x x \V) = {l + 2K 1 x 1 ) v Y! 

= (l + 2K lXl ) v E{z e ,z fl V), (6) 




FIG. 1: A configuration of chains placed on the lattice. Inter- 
nal monomers are represented as white circles and endpoint 
monomers are black circles. The statistical weight of this con- 
figuration with 5 chains is z\ z]P . 



tice has V sites, the density of endpoint monomers is p e = 
N e /V, the density of internal monomers is p t = Ni/V, 
and the total density of monomers is p = p e + Pi- The 
densities may be obtained from the partition function: 



and 



Pe 



E dzi ' 



z e 

E dz e 



(9) 



(10) 



Let us now consider the model defined on a Cayley 
tree with coordination number q. Due to the hierarchical 
structure of this lattice, we may define partial partition 
functions on rooted subtrees, and it is not difficult to 
write down recursion relations for subtrees with an ad- 
ditional generation [13j |. For the present model, we may 
define two partial partition functions go and g%, with- 
out and with a polymer bond on the root bond of the 
subtree, respectively. Considering the operation of at- 
taching a = q — 1 n-generations subtrees to a new site 
and bond, we obtain recursion relations for the partial 
partition functions of a subtree with n + 1 generations. 
The contributions to this recursion relations are shown 



3 



in Fig. [2 and the resulting expressions are 



and 



9a = 9o + az e9 9i + 



9i = z c.go + vZi9o 1 9i- (12) 
If we now define the ratio of partial partition functions 



Root bond 



FIG. 2: Contributions to the recursion relations for the partial 
partition functions. Fhe first line are the contributions to g' Q 
and the second line shows the contributions to g[. 

R = gi/go-, we will find the following recursion relation: 

z e + aziR 



R' 



1 + <TZ e R 



g(g-l) 



z l R 2 



(13) 



In the thermodynamic limit, the fixed point of this re- 
cursion relation is reached. 

The partition function of the model on the Cayley tree 
may be obtained if we consider the operation of attaching 
q subtrees to the central site of the tree. The result is: 

" 9 i 9-1 i 1\1 ~ 1) <?- 2 2 ft a\ 

- = 9a+qze9l 3i + z i9 l g v (14) 

The Bcthc lattice corresponds to the central region of the 
Cayley tree, and therefore we concentrate our attention 
on the densities of endpoint monomers and of internal 
monomers on this site. The results are: 



qz e R 



1 + qz e R ■ 



g(g-i) 



z t R 2 



and 



g(q-i) 



ft 



z t R 2 



l + qz e R+siSzll Zi R2 



(15) 



(16) 



For given values of the activities z e and Zi, the fixed point 
value of the ratio R may be found using Eq. I[13p and 
thus the corresponding densities of endpoint and internal 
monomers at the central site of the tree are defined by 
Eqs. (fT5|) and (jTHJ), respectively. For z e > 0, the densities 
are smooth functions of the activities, but in the polymer 
limit (infinite chain) z e — > a continuous phase transition 
is found between a non-polymerized phase ft = and a 
polymerized phase ft > 0, at Zi — l/(q — I) fhH. [l5j. 

If we call s(p e ,pi) the entropy per site of the model, 
we notice that the ratios of the chemical potentials and 
fcgT are given by 



Me 



= In z e = - —— 



ds 
dp. 



(17) 



In Zi 



ds_ 

dpi 



(18) 



We then may obtain the entropy integrating the chemical 
potentials. To do this, we start inverting equations for 
the densities (IT51) and ([TBI), which leads to: 



and 



q(l-pi-p e )R 

2^ 

q(q-l)(l- Pl -p e )R 2 



(19) 



(20) 



If these activities are now substituted into the fixed point 
equation R' — R, we obtain, from Eq. (JT3J) the fixed 
point value of the ratio: 



R 



2ft 



q-2pi 



(21) 



Substituting this result in Eqs. (fTTj)) and (|2"0")) . we obtain 
the activities as functions of the densities in the thermo- 
dynamic limit: 



Pe.\Jq - 2pi — Pe 



and 



q(l - Pi- p e )V 2 Pi + Pe 

2p t (g - 2p, - p e ) 

9(9 - 1)(1 - Pi - Pe){2Pi + Pe) ' 



(22) 



(23) 



The entropy may now be calculated integrating the Eqs. 
(fT7)l and ([T8j) in the (p e ,pi) plane. We choose the trajec- 
tory (0,0) -> (0,pi) -> (p e ,Pi), so that: 



s(p e ,Pi 



In 



q~2p 



g(?-i)(i-p) 



dp — 



hi 



pyjq - 2pi - p 



q(l - pi - p)^2p l + p 



dp. (24) 



The result of this integration is: 

s(p e ,pi) = --(q-2pi)]nq + pi]n(q- 1) - p e ln — 
2 q 

(q - 2pi - p e ) \n(q - 2p { - p e ) - 



1 



(1 - p t - p e ) ln(I - Pi - p e ) + 

-(2ft + pe) ln(2ft + p e ) - ft m(2pi). (25) 

The entropy may be compared to results in the liter- 
ature in some particular cases. For p e — the case of 
infinite polymers is obtained, and the expression reduces 
to Eq. (26) in reference [5]. The entropy of dimers on 
the Bethe lattice (Eq. (27) in reference @, where in the 
correct result pz should be replaced by p/z) is obtained 
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for pi = 0. Finally, the one-dimensional result (Eq. (9) 
in reference [l2|) is recovered for q = 2. To compare 
the result with the monodisperse case, it is convenient to 
rewrite the entropy Eq. (|25[) as a function of the mean 
molecular weight M = 2(1 + pi/ p e ) and the total density 
of monomers p = pi + p e . The resulting expression is: 



Sm(p) 



p M — 



2 1 ( i\ 2 Pi 2 P 



q m-i \ 



2p- 



M- 1 
q~M 



(l-p)]n(l-p)+p 



M 



M 



1 , ( ~ M 
-In 2p 



M 



M - 



P- 



M 



M 



(26) 



We may now find the contribution to the entropy of poly- 
dispersivity, calculating Asm(p) = Sm(p) — s m(p) for 
M = M, where sm{p) is the entropy of the monodis- 
perse case (Eq. (22) of reference Q). We find a result 
which is independent of q: 

As M (p) = -^[(M - 1) ln(M - 1) - (M — 2) ln(M - 2)], 

(27) 

and therefore this result is identical to the one obtained 
in the one-dimensional case studied in reference fl2| . 

Finally, we may find the distribution of the molecular 
weights in the polydisperse case. To obtain this result, we 
must keep track of the number of monomers incorporated 
into each chain, and we may define multiple subtree par- 
tial partition functions gM, M = 0,1, . . ., in such a way 
that go corresponds to a subtree with no polymer bond 
at the root, as before, and gM stands for a subtree with 
a polymer bond on the root connected to M monomers 
above. The recursion relations for these partial partition 
functions will be: 



So = 9o + az e9o 1 2^ 9M + 

M=l 

2 Zi 9o 2^ 9M9N, (28) 

M,N=1 



fo- 1 E 



9l = z e9 , 



9m = ° z i9o 9m-i,M = 2,3,. 



(29) 
(30) 



We proceed defining the ratios of partial partition func- 
tions Rm = ghi/go, for M — 1,2,.... The recursion 
relations for the ratios will be: 



R' M - °^^,M = 2,3, 



D 



(31) 
(32) 



where 



D = 1 + az e Rm + a[a ~ ' zi r mRn- (33) 



Inspection of these set of recursion relations leads to the 
following Ansatz for the fixed point values of the ratios 
Rm — a M ~ l R\, where a — R\azijz e . The ratios Rm 
are related to the ratio R defined above through R = 

Em=i r m = Ri Em=i " m_1 = Ri/i 1 ~ «)■ Usin S the 

fixed point value for R found in Eq. (|2ip . we may then 
find the values of R\ and a in the thermodynamic limit. 
The results are: 



Ri 



Pe 



yj{p e + 2pi)(q- p e - 2pi)' 



and 



2 Pl 



Pe + 2p t 



(34) 



(35) 



Now we consider the operation of attaching q subtrees 
to the central site of the tree. The probability that an 
endpoint of a chain with M monomers is located on the 
central site will be equal to qz e RM-i, and the probability 
to have an endpoint of a chain located on the central site 
is equal to qz e R. Therefore, the probability to find a 
chain with exactly M monomers among all chains on the 
lattice with mean molecular weight M will be given by: 



R 



M-i 
~R~ 
1 



(1 — a)a 



M-2 



M — 1 \M —1 



M 



M-2 



(36) 



This exponential distribution of molecular weights is 
again independent of the coordination number q and 
therefore is identical to the one obtained in [12j for the 
one-dimensional case q — 2. 



III. FINAL DISCUSSIONS AND COMMENTS 

Although the equilibrium polymerization model was 
extensively studied in the literature, the main focus of 
these studies was its critical behavior, which occurs in 
the limit K\ — > (z e — > 0), where the mean molecular 
weight of the chains diverges. Here, similarly to what was 
done by Dudowicz et al [9( for a similar model, which es- 
sential differs from the present one only by the fact that 
it is canonical with respect to the iniciator molecules, 
we turn our attention to the region of parameters where 
the model is non critical. The calculations performed 
by Dudowicz et al were in the Flory-Huggins approxima- 
tion, which for simple polymer models is equivalent to the 
Bethe lattice solution. It is therefore not surprising that 
boths calculations lead to an exponential distribution of 
molecular weights of the chains. The model studied by 
Dudowicz et al is canononical with respect to the initia- 
tor molecules, while the present model is grand-canonical 
with respect to the endpoint monomers. In the polymer 
limit, the distribution of molecular weights should in gen- 
eral be described by the scaling behavior [16| 



M=l 



M,N=1 



t m ~ M 7_1 exp(-M), 



(37) 
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in the limit of small overlap, reducing to an exponential 
decay at large overlaps. The critical exponent 7 is in gen- 
eral larger than 1, but its classical value is 1. Since the 
Flory-Huggins approximation and Bethe lattice calcula- 
tions should lead to classical exponents, the exponential 
decay of the distribution of molecular weights is expected 
in the polymer limit M — > 00, but nonexponential decay 
would be allowed for finite M. We are presently inves- 
tigating if this actually occurs when the model is studied 
using better (mean field) approximations. 

Although the details of the calculations above were 
presented for the restricted equilibrium polymerization 
model, without the presence of one-site chains, it is easy 
to obtain the correpondin results for the original model, 
using Eqs. |0 and |[8|. For example, the density of 
sites occupied by the one-site polymers is 

Pl= Yd^> (38) 
making xi — 1/2. The result is 

Pl = T ^(l-Pe- P i). (39) 



A similar calculation may be done for the density of poly- 
mers (including the one-site chains) p p = N p /N, leading 
to: 

*=i+k( 1+ ^*r*-"\) (40) 

while the density of bonds is given by: 

Pb = y + Pi- (41) 
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